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Under certain circumstances, a three-body decay width can be approximated by an integral involving 
a product of two off-shell two-body decay widths. This "angle-average" approximation is used to 
calculate the tyNN decay width of the d*(J^ = 3 + ,T = 0) dibaryon in a simple A 2 model for the 
most important Feynman diagrams describing pion emissions with baryon-baryon recoil and meson 
retardation. The decay width is found to be about 0.006 (0.07, 0.5) MeV at the d* mass of 2065 
(2100, 2150) MeV for input dynamics derived from the Full Bonn potential. The smallness of this 
width is qualitatively understood as the result of the three-body decay being "third forbidden". 
The concept of i forbiddenness and the threshold behavior of a three-body decay are further studied 
in connection with the irNN decay of the dibaryon d 1 (J 7 * = 0~,T = 0or2) where the idea of 
unfavorness has to be introduced. The implications of these results are briefly discussed. 



PACS number(s): 

I. INTRODUCTION 

Almost all theoretical dibaryons with exotic quark 
structures have masses above the irNN threshold |Q. 
They can thus decay into irNN states. From the view- 
point of dibaryon searches, the most promising candi- 
dates should have the narrowest widths, otherwise they 
do not stand out clearly above the background. However, 
their irNN widths should not be too small, otherwise too 
few pions will be available to help with the identification. 
Hence a qualitative understanding of the irNN decays of 
these dibaryons is of considerable interest when contem- 
plating an experimental search for these dibaryons. 

A particularly interesting dibaryon is the d* , of quan- 
tum numbers J^T = 3 + 0. The interest comes from the 
possibility that its mass might be unusually low, thereby 
indicating an unusual structure or dynamics The 
purpose of this paper is to give rough estimates of its 
irNN decay width when it is treated approximately as 
an S-wave bound didelta state of intrinsic spin S = 3. 

The irNN decay of this model d* cannot occur via 
pion emission from one of the constituent A's because 
the spectator baryon will remain a A. Hence the baryons 
must interact at least once to turn the spectator into a 
final-state nucleon. This is true too for pion emission 
from a virtual meson which must be made to appear in 
the system. With one pion emission vertex and two in- 
teraction vertices for the spectator interaction (referred 
to below as a recoil), the leading Feynman diagrams are 
processes containing three vertices, as shown in Figs. 1 
and 2. 

Of the processes shown, Fig. la-c are recoil diagrams 
describing ir emission by interacting baryons. Fig. ld- 
f are retardation diagrams describing emission from a 
baryon when a virtual meson is in the air. Figs, lg-i 
are some of the relativistic corrections coming from NN 
pairs. Finally Fig. 2 gives the contributions for emission 
from an exchanged meson itself. We shall calculate the 
decay width for the most important of these processes 



with the help of an approximate angle-average formula 
developed in this paper. 

This paper is organized as follows: The general na- 
ture of the irNN decay of the model d* is described in 
Sec. II for P- and S-wave pion emissions. The concept 
of ^-forbiddenness for describing the threshold behavior 
of three-body decays is introduced in analogy to that for 
nuclear (3 decays. Both P-wave and S-wave pion emis- 
sions from baryons described by Fig. 1 are found to be 
"third-forbidden" . The processes shown in Fig. 2 for 
pion emission from virtual mesons in the leading order 
are found to vanish for both P-wave and S-wave pion 
emissions. 

In Sec. Ill, an angle-average approximation for the 
width of a three-body decay of the type shown in Fig. 
1 is obtained in the form of a sum of integrals over the 
product of two sets of off-shell two-body decay ampli- 
tudes. Numerical results are given in Sec. IV where this 
angle-average approximation is applied to Figs, la-c, us- 
ing baryon dynamics derived from the Full Bonn poten- 
tial |f| and from on-shell NN i-matrices constructed from 
experimental NN phase shifts ||. In Sec. V, retardation 
and pair contributions are included as well. 

Sec. VI gives a discussion of the threshold behaviors 
of irNN decays, which are controlled but not solely de- 
termined by the £-forbiddenness of the decay. The decay 
of the d'(,r = Q~,T = 0or2) dibaryon is studied to il- 
lustrate the complications caused by decays with abnor- 
mally large power dependences on the external momenta 
near threshold. 

Finally, brief concluding remarks are made in Sec. 
VII, where the implications of our results for the quark- 
delocalization model of d* |2]— Q are touched upon. In 
particular, the calculated decay width of 70 keV at 
to* = 2100 MeV seems to imply that there might be 
too few decay pions at low dibaryon masses to help with 
particle identification. 
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II. GENERAL CONSIDERATIONS 

The dibaryon d* is an isoscalar, high-spin (J 71 " = 3 + ) 
state. When treated as a A 2 bound state, its dominant 
spatial component is 7 S3. Fig. 1 gives the lowest-order 
Feynman diagrams for its decay into the ttNN final state 
by pion emission from a baryon. For P-wave pion emis- 
sion, the NN system in the ttNN final state has ST = 01, 
the possible partial waves being 1 D 2 and 1 G&. Hence the 
7riViV decay of d* through this dominant 7 S3 component 
requires three units of baryon spin change, i.e. a rank-3 
intrinsic spin tensor. 

This rank-3 spin tensor must be combined into a rota- 
tional scalar as the perturbation responsible for the de- 
cay by scalar multiplication with a rank-3 spatial tensor 
made up of internal or external momenta. The Feynman 
diagrams shown in Fig. 1 describe some of the leading- 
order pion emission by a baryon in our simple model of 
d* . The dominant process involves a P-wave pion emis- 
sion that depends linearly on the momentum P3 = of 
the emitted pion. This means that the remaining qq or 
BB interaction must generate a rank-two spatial tensor 
to trigger the decay. This cannot be done with a rank-0 
central interaction, or a rank-1 spin-orbit interaction. A 
rank-2 tensor-force interaction is needed. Furthermore, 
the final NN state with the lowest orbital angular mo- 
mentum that appears from the initial dominant spatial 
7 S 3 -state of d* is 1 D 2 . 

To the extent that d* is dominated by the A 2 com- 
ponent containing two t = 3/2 constituents, the tensor 
interaction for recoil must convert at least one of these 
constituents into a ( = 1/2 nucleon. Thus the tensor 
force must be an isovector interaction that arises, for ex- 
ample, from the exchange of an isovector meson such as 
7r or p. This means that the operator in the BB inter- 
action responsible for the decay in our simple A 2 model 
of d* is (<Tj x (Tj)W(Ti.Tj), where i and j refer to two 
separate baryons, or to quarks in these baryons. The op- 
erators <r and r are Pauli spin operators for quarks, and 
generalized spin operators for baryons . 

It is useful to enumerate the possible components of 
both the initial d* and the B 2 states appearing after P- 
or S-wave pion emission. We shall use a baryon nota- 
tion with baryon components restricted to only A and 
N because these baryons have the lowest masses and the 
same spatial s 3 quark structure. They can therefore be 
expected to mix strongly with one another. Under the 
circumstances, the available states are 

d*{.r = 3+,T = 0) : A 2 [ 7 {S,D,G,I)}; 

N 2 [ 3 (D,G)}. 

B 2 (.r = 2+/4+,T= 1) : N 2 ^(D,G)]; 

AN[ 5 (S,D,G,I)],[ 3 (D,G)]; 
A 2 [ 5 (S,D,G,I)], [ 1 (D,G)]. 
B 2 (,r =3~,T= 1) : N 2 [ 3 F}; 



AN[ 5 (P,F,H)],[ 3 (P,F)}; 
A 2 [ 7 (-P, F, H)}, [ 3 (P, F)]. (1) 

The decay processes shown in Figs, la-c can then be de- 
scribed by certain decay equations. For P-wave emission, 
they are 

A 2 ( 7 S 3 ) - N 2 ( 3 D 3 ) -> N 2 CD 2 ) + tt(4 = 1); 
A 2 ( 7 S 3 )^ AN( 5 S 2 )+n(£^ = l), 

AN( 5 S 2 ) -> N 2 CD 2 ); 
A 2 ( 7 S 3 )-> A 2 ( 5 S 2 )+n(£„ = l), 

A 2 ( 5 5 2 ) -» N 2 CD 2 ) . 



(2) 



In this paper, we are particularly interested in low- 
mass dibaryons decaying close to its ttNN threshold. 
Such near-threshold decays show certain general kine- 
matical features reminiscent of nuclear (3 decays, features 
that originate from the dominance of centrifugal poten- 
tial barriers on the outgoing decay products, and are con- 
trolled by their orbital angular momenta. The transition 
amplitude can be expanded in a Taylor series in these 
final-state momenta of which only two, say of the 
pion and pjv of one of the nucleons, are independent. If 
£ w is the pion orbital angular momentum in the dibaryon 
cm. frame, and £m is the relative orbital angular mo- 
mentum of the TV TV system in the final state, the leading 
terms in the Taylor expansion are necessarily of the type 
3 ; ^m^(P7r)3 7 £ Jv m J v(PAr), where 3^ m (p) is a solid spheri- 
cal harmonic. The decay width near threshold is then 
proportional to (p l Jp 1 ^) 2 together with additional fac- 
tors coming from the three-body phase space which we 



shall find in Sec. |VI| to be p 2 p%- We shall call this decay 
I- forbidden when £ — £ n + £n is nonzero. With £ n = 1 
and £n — 2, the P-wave pion decays of Eq. (Jf2]) are all 
third- forbidden. 

This concept of £-forbiddenness is familiar not only 
from nuclear decays, but also from the simpler pro- 
cesses of two-body decays. For example, the decay of the 
A resonance involves P-wave pion emission, of £ n = 1. 
Hence the decay is £ = 1 forbidden in our language. The 
two-body phase space gives an additional factor of p v . 
The decay width is then roughly proportional to p 3 . This 
is a well-known result. 

The ^-forbiddenness of the decay d* — > irNN is de- 
termined as follows: The pion has both orbital and total 
angular momenta and parity 7iv = (— The ini- 
tial state has quantum numbers Jf — 3 + , T — 0. Hence 
the quantum numbers of the most favorable final NN 
state is T = 1, £n = \ J% — £^\, and irpj — WiTT T . From 
these selection rules, we can determine that the decay 
d* -> irNN is £ = Ji = 3 forbidden for £„ = - 3. For 
higher £„ > Ji, it is £ = 2£ v — Ji forbidden. 

These selection rules for ttNN decays have been spec- 
ified in terms of the relative coordinate of the final NN 
state, i.e. a Jacobi coordinate. This is the most natu- 
ral choice because wave-function antisymmetrization and 
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final-state interactions can then be described most con- 
veniently. However, it is not the only one possible. Se- 
lection rules can also be specified in terms of the one 
of the final-state nucleon coordinates. This choice could 
occasionally be more convenient, for example when the 
nucleon is a real spectator. An example of this choice 



will be given in Sec. VI 



Let us now return to the decay of d* . For S-wave pion 
emissions, the final A A system should have odd spatial 
parity, and therefore ST =11. Only one A 2 partial wave 
is then possible, namely 3 F 3 . Figs, la-c for S-wave pion 
emission from a baryon are described by the decays 

A 2 ( 7 S 3 ) - A 2 ( 3 P 3 ) -» N 2 ( 3 F 3 ) + tt(4 = 0); 
A 2 ( 7 S 3 )-+ AA( 5 P 3 )+^(C =0), 

AA( 5 P 3 ) -» N 2 ( 3 F 3 ), 
A 2 ( 7 S 3 )^ A 2 ( 7 P 3 )+^ = 0), 

A 2 ( 7 P 3 ) -> A 2 ( 3 P 3 ) , (3) 

respectively. Like P-wave pion-emission amplitudes, the 
S-wave amplitudes interfere among themselves, but they 
add only incoherently to the P-wave contributions to the 
decay width. 

Now the basic S-wave pion emission vertex is known 
H to be of order w/m relative to P-wave pion emission, 
where the baryon (or quark) energy transfer w — E — E' 
is expected to be significantly smaller than the baryon 
mass m for low-mass d* . Given the fact that these S-wave 
decays are also I — 3 forbidden, we conclude that they 
can be neglected in comparison with P-wave emissions, 
at least for the rough estimates attempted in this paper. 

Other types of S-wave pion emission are possible if the 
Feynman diagrams involve orbitally or radially excited 
baryons A* or A* as well. These additional baryon con- 
figurations can be expected to be even less important 
because the mBB* coupling constants are much smaller 
than the meson coupling constants to A or A because of 
the spatial quark excitations. 

Among higher-order diagrams not included in Fig. 1 
are those decribing final-state interactions between the 
outgoing nucleons. They are not expected to be very im- 
portant because the nucleons have large relative orbital 
angular momenta in our final states. Thus re-scattering 
effects in d* — > ttNN decays could be quite different from 
those of threshold pion production in nuclear reactions 

Let us turn next to Fig. 2 describing pion emission 
from virtual 7r — p mesons. The operators involved in Fig. 
2a have been derived in connection with the study of inrp- 
exchange three-nucleon forces JTT| . Two pion-emission 
operators appear: The first is (er 2 .q)(T2 x ti).$, where 
q is the momentum of the virtual pion and <fr is its 
wave function. The spin-isospin operators appearing here 
could be those of quarks or of baryons (including tran- 
sition spins) depending on the dynamical model used. 



For simplicity, however, we shall visualize them in this 
section as generalized baryon spin operators. 

This pion-emission operator describes S-wave pion 
emission. For the two diagrams in Fig. 2 together, 
S-wave pion emission involves the combination (<Ti + 
<T2).q(xi x T2).3». This is a rank-1 spin operator. Hence 
the direct decay of the major component A 2 ( 7 5' 3 ) of d* 
(i.e. Fig. 2) to the only available A 2 final state of 3 F 3 
is not possible. The process has to go through a minor 
component of cither d* or the final B 2 state. A minor 
component of the final B 2 state is one which contains one 
or more A's that must be converted to nucleons before 
the decay is completed. 

The second operator appearing in Fig. 2a is propor- 
tional to (<T2.q)c 1 .(p 3 x q)(i~2 x t%).&, obtainable also 
from 7T7T/9-exchange three-nucleon forces. It gives rise to 
a pion emission operator of the form — 2(<ti x <x 2 ).[q x 
(P3 x l)](' r i x T 2)-&, describing P-wave pion emission. 
The intrinsic spin operator appearing in it is also only 
rank 1. It too cannot connect the dominant A 2 ( 7 5 3 ) 
component of d* to the available A 2 final states 1 £>2 or 
1 G*4. Hence P-wave pion emision by a virtual meson is 
again possible only through a minor component in the 
initial or the final state. 

To summarize, to the lowest order in the strong- 
interaction vertices, the decay d* — > irNN involves the 
processes described by Fig. 1 involving pion emissions 
from baryons. 



III. ANGLE-AVERAGE APPROXIMATION FOR 
THREE-BODY DECAY WIDTHS 

Estimates of the width of a three-body decay of the 
type shown in Fig. 1 are greatly facilitated by using an 
angle-average approximation developed in this section. 

We begin by noting that Fig. la can be interpreted as 
the off-shell decay A — > irN of a small NN( 3 D 3 ) compo- 
nent of d* that has one spectator nucleon on its energy 
shell. This component is generated perturbatively by the 
recoil interaction from the A 2 ( 7 S 3 ) state. For this reason, 
the energy denominator that appears is the bound-state 
expression 



D a = m* - 2P(p 2 ), 



(4) 



where m* is the d* mass and E is a nucleon energy in 
the d* cm. frame. The f-matrix appearing in the recoil 
interaction is an off-shell i-matrix between two baryons 
calculated at the initial-state energy i/s = m*. 

On the other hand, Figs, lb and c involve pion emis- 
sion into abnormal or minor A 2 or AA components of 
the final B 2 channels. These abnormal components are 
off-shell with finite spatial extensions when the available 
energy is below their breakup threshold. 

Figs, lb and c are also similar in that their energy 
denominators are scattering-state quantities of the type 
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D b = E( Pl ) + E(p 2 ) - E A (-p) - E(p - p 3 ), (5) 

D c = E{ Pl ) + E(p 2 ) - E A (-p) -E A (p- p 3 ), (6) 

where E A is a A energy in the d* cm. frame. Nei- 
ther energy denominator vanishes below the appropriate 
threshold. The off-shell BB f-matrix is that calculated 
at the energy of the final-state NN pair. 

The n-body decay of an object of mass m* can be 
expressed conveniently by the Fermi golden rule [pL 2|| : 



r = (2tt) J 



1 

dE 



n 

i=2 



d 3 P 



-fffi| 2 )spin 



(2tt) 3 
x5(m* - E)dE, 



(7) 



where H' is the effective perturbing Hamiltonian. The 
n-body density of states that appears can be expressed 
nominally in terms of a product of densities of states of 
fewer bodies. If the spin-average squared matrix element 
(|-ff^| 2 ) S pin is also factorable into appropriate factors re- 
ferring to fewer bodies in the problem, the integral could 
be simplified to make the physics more transparent. 

To study this possibility for d*, we first note that 
each Feynman diagram of Fig. 1 has the simple form of 
H 2 D~ 1 H^, where Hi is one of the interactions and D is 
the energy denominator between them. The three-body 
phase space is unfortunately made complicated by the 
requirement of energy conservation which leaves behind 
the integrations over 



^ ^ d3 P* = E^E^dE^n^d^ . (8) 



Here use has been made of the momentum-conserving 
relation 



E\ = p\ + p\ + 2p 2 P3 cos 6> 23 



in the form 



d COS #23 Ei 



dEi 



P2P3 



(9) 



(10) 



Thus only one set of angle integrations, here d 2 f2 3 , is 
formally identical to that for the decay in free space. 

What is left of the second set of angle integrations, 
namely d</>23, has been partly changed by the requirement 
of energy-momentum conservation for all three particles. 
As a result, kinematical quantities such as E\ and 023 are 
no longer independent of each other. Many of the angu- 
lar functions in the integrand become quite complicated. 
The energy denominator D too could be a function of 
some of these angle variables. 

We are interested here in an angle-average approxi- 
mation obtained by adding an extra angle integration 
(l/2)d cos 6*23 to Eq. (||) to restore the full angle integra- 
tion of d 2 r^2- The idea is to undo the angular correlation 



between the matrix elements of the interactions Hi and 
H 2 (so that the the angle integrations d 2 pi can be done 
independently), but keep the proper three-body phase 
space of the Dalitz plot |l^,[l3| . 

Before this can be done, another complication has to 
be taken care of. The integrand (|(^f')fl| 2 )spin contains 
the propagator factor D~ 2 which depends on the external 
hadron energies Ei and certain internal baryon energies. 
The external hadron energies are angle-independent, but 
some of the internal baryon energies do depend on the 
integration angles. If such angle-dependent internal en- 
ergies are replaced by suitable angle-averaged values, for 
example by the replacement 



P3 



pI, 



(ii) 



the factor D~ 2 simplifies to an angle-independent ap- 
proximant (D~ 2 ) that can be taken out of the angle in- 
tegrations. 

The remaining factor in the integrand is now 



(i( J ff 2 i/3) f ii 2 > S pi„=n(i(^) fi i 2 ) s P i - 



(12) 



where the state labels i, f in are the initial- and 

final-state labels appropriate for the interaction Hi . The 
two sets of angle integrations can now be done indepen- 
dently, leading to the "off-shell" two-body decay widths 



Ti(Wi] 



pin, angle ; (13) 



{2ir) 2 dWi 



where 



(|(FOfir) S pin,anglc = ^ J ( | (#i )fl | 2 ) spin^Pi • (14) 

These decay widths are in general off-shell because the 
energy in each final state differs in general from that in 
the initital state. The effective energy Wi in the density 
of states remains to be chosen. 

Eq. (M) has thus been simplified to 



-!- / dW 2 dW z 

Z7T 



Ex 



%P2P3 , 

r2(W2)CD- 2 )r 3 (W3). 



(15) 



This is the basic angle-average approximation used in this 
paper. In this expression, the effective energies Wi are ar- 
bitrary integration variables to be chosen for convenience 
since the integral is actually independent of their choices. 
Whatever the choice, the remaining integrations should 
be performed over an appropriate Dalitz area that repro- 
duces what is left of the three-body phase space after the 
angle integrations. 

We choose W3 = E$ so that T3 is actually the static 
limit [where <ti m(nucleon)] of the physical width 
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for pion emission. We use W 2 = 2E 2 in order to make 
the width T 2 for Fig. la an off-shell d* — > NN decay 
width with a two-body density of states defined for the 
iV 2 system of total energy 2E 2 . 

The actual nNN decay amplitude is of course a 
sum over contributions from several Feynman diagrams, 
though limited in this section for convenience in presen- 
tation to only Figs. la-c. To treat their interference, we 
shall need the complex transition amplitudes for the BB 
recoil 



where 



F, 



(16) 



The final decay width will then be the sum of contribu- 
tions 

r = y"T(a, /?), with a, (3 = a, b, ore , (17) 

a,0 



where 



T(a,/3) 



•in 



P2P3J 
xF£ a F 2l3 {(D a D )- 



(18) 



if there were no additional complications in the angle 
and spin averages for the "off-diagonal" terms. An extra 
factor of 4 appears because there are two A's in d* . The 
diagonal terms T(a, a) are positive, while the interference 
terms T(a, (3 ^ a) could be negative. 

Further development of this formula is possible under 
special circumstances. Matrix elements at corresponding 
vertices in different diagrams can be related to one an- 
other through the quark model [p]Jl4|. For example, the 
irqq vertex for pion emission from quark i can be taken 
to be the NR expression 



f-rrqq 



(CT;.P3)(T;.*) , 



(19) 



where $ is the pion wave function, and f„ qq = (3/5)/ w atat 
is the coupling constant. The NR quark model can then 
be used to relate the B — > irN decay widths that appear 
in these Born diagrams to just the basic width 



r 3b (p 3 ) = r(A -> Nir) 



(20) 



together with additional factors that represent changes 
in the reduced matrix elements of operators appearing in 
the interaction. Additional details will be given in the 
Appendix. Re-scattering effects for the emitted pion will 
not be included in the present study. 

The BB — ► NN transition is more complicated. The 
transition matrix element (H2) fi in Figs, lb and c should 
be calculated in the d* rest frame, not in the final-state 
N 2 cm. frame where the dynamics is usually specified. 
This means that 



F 2l (W 2 ) = 5 p T 2l {p* 2 ), 



(21) 




(22) 



and T is an off-shell transition amplitude containing the 
two-body density of states calculated in the N 2 cm. 
frame where the nucleon momentum is 



P2 = \l J m i2 



mr 



and 



m 12 = (E 1 + E 2 f 



(23) 



(24) 



is the invariant mass of the A^ 2 system. 

We could still use the NR quark model to relate the 
different operator matrix elements of different diagrams 
so that the same basic transition amplitude T 2a appears 
in different diagrams, albeit evaluated at different ener- 
gies and with different additional factors. Finally, spin 
averages must be performed term by term for Eq. JTt|), 
as discussed in the Appendix. In this way, we get the 
final result for Figs, la-c of 



r «il 

aC 127T 



dE 2 dE 3 



D„ 



P2P3J 
2T 2a {p* 2 ) 
3D b 



^2a{P* 2 ) 

3D C 



(25) 



involving angle-averaged quantities. This is the decay- 
width expression actually used in our calculations. 

We could have used empirical mBB coupling constants 
or some other model of mBB dynamics (l4| instead 
of the NR quark model to calculate different BB recoil 
amplitudes. The result will differ somewhat from those 
shown in the equation. However, the differences are not 
expected to be very large from the viewpoint of the rough 
estimates attempted here. So we shall not consider these 
alternatives in this paper. 

Equation ( p5| ) has a particularly simple structure if the 
recoil amplitude T 2a is calculated in the Born approxi- 
mation or evaluated at a common nucleon momentum p 2 
instead of at two different values. On ignoring the small 
difference in the density of states in different N 2 frames 
as well, wc find 



^lmdE s [^-) 



1 



+ 



P2V3J 
2 1 



r3fcF2a 
2 



where 



2a 



D a 3D b 3D C 



r(d* -> a 2 ( 1 d 2 )) 



(26) 



(27) 
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This very rough formula for T shows that the relative 
importance of the diagrams is approximately controlled 
by the energy denominators Di. 

To estimate these energy denominators, we note that 
an external nucleon has the median energy of 



E(p 2 ) « Um N +E n 



961 (988) MeV, 



(28) 



where the numerical value is for m* = 2100 (2200) MeV. 
For small d* masses, the internal baryon energies that 
appear are dominated by the internal momentum p of d* 
of radius r* (usually taken to be 0.7 fm): 



<p 2 > 



9 



16r* 



(210MeV/c) 2 , 



N 



-< p 2> = 962 MeV, 



E(p) 



E A (p) « v/ m A+<P 2 > = 1250 MeV, 



(29) 



(30) 



(31) 



and Ei(p — p 3 ) w Ei(p). These gives the rough estimates 
at m* = 2100 (2200) MeV of 

D a w 178 (224) MeV; 
D b w -290 (-236) MeV, D c w -578 (-524) MeV. (32) 

These approximate energy denominators can be used 
in Eq. (^) to show that the decay widths from the three 
diagrams taken individually are roughly in the ratio 

r(o, a) : r(6, b) : T(c, c) w 1 : (-0.41) 2 : (-0.10) 2 (33) 

at 2100 MeV, and 

T(a, a) : T(b, b) : T(c, c) » 1 : (-0.63) 2 : (-0.14) 2 (34) 

at 2200 MeV. Because the denominator of Fig. la has 
sign opposite to those in Figs, lb and c, the amplitude 
from this dominant process interferes destructively with 
those from the latter diagrams, as shown explicitly in the 
last two equations. This destructive interference leads to 
a total decay width for Figs, la-c that is roughly an order 
of magnitude smaller than T a from Fig. la alone. 



IV. RESULTS FOR D* 



■kNN 



To obtain a wNN decay width T, we need to integrate 
Eq. (^5|) over the three-body phase space left over after 
angle integrations. It seems useful to give a rough order- 
of-magnitude estimate of the dominant term r a from Fig. 
la as we go over some of the technical details involved in 
the integration. 

The two-body width for P-wave pion emission is best 
expressed in terms of the experimental decay width Ta = 



T(A -> nN) of 120 MeV for a free A of mass m A = 1232 
MeV: 



r 3b = r A (|||V^ 2 ^-,i) 



m 2 ^ 



(35) 



where the starred parameters are in the A cm. frame: 



227 MeV, E 



N 



966 MeV 



(36) 



The Gaussian factor comes from a baryon formfactor cal- 
culated by using Gaussian wave functions in a NR quark 
model. It depends on a parameter 



a = r 2 /6 = 1.54 x 10" 6 MeV" 



(37) 



related to the proton radius r p . The value of r p = 0.6 
fm used here is the one commonly used in many NR 
quark models, including It is not the experimental 
proton charge radius. Note that the pion momentum p^ 
appearing in Eq.(|35"|) is that in the d* rest frame, and 
that the distortion of the outgoing pion wave function 
has been neglected. 

Momentarily ignoring the Gaussian factor from the 
baryon formfactor, we find for m* = 2100 (2200) MeV 



L 3b max — 

59 (256) MeV, 

\ _ r3fr max 

P3 



2.5p. 



3 max 



0.14(0.38) ; 



(38) 



at the maximum value P3 ma x = 165 (270) MeV of p%. 
The average shown in the last expression involving a de- 
nominator factor of 2.5 is obtained as follows: For non- 
relativistic (extremely relativistic) kinematics, the phase 
space is roughly rectangular in p\ (P3). The average then 
has a denominator factor of 2 (3). For most of the pion 
energies studied here, the kinematics is neither nonrela- 
tivistic and extremely relativistic. So we take a denomi- 
nator factor of 2.5. 

The two-body width T2a(p2) for the two-body decay 
d* — » NN is much less familiar, but it has recently been 
studied in Q. It is given by Eq. (38) of that paper, an 
expression that works for any final nucleon momentum 
P2 in the rest frame of the final-state N 2 system. It can 
be used for the off-shell momenta appearing in our dia- 
grams, provided that the reduced mass n*j appearing in 

the equation is taken to be E%/2 (where E2 — \/p\ + m2 ) 
and not m*/4. 

Different dynamical inputs are possible. We use the 
full Bonn (FB) potential in the Born approximation, a 
choice referred to below as the FB model. Of course, the 
Born approximation of a static potential is rather crude 
because it contains no off-shell effect, but it represents 
a familiar starting point. The resulting width parameter 
^2a(p) /p 2 is shown by a solid curve in Fig. 3 as a function 
of the nucleon momentum p in the cm. frame of the final 
NN state. 
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Certain re-scattering effects not included in Fig. 1 
could be included by using empirical NN t-matrices con- 
structed by Love and Franey (LF) || from experimen- 
tal NN phase shifts. These empirical t-matrices are of 
course on-shcll NN t-matrices, but here evaluated at 
the same momentum transfer q as the needed off-shell 
t-matrices. So no dynamical off-shell effects can be in- 
cluded here. To my best knowledge, there is no sim- 
ple way to extrapolate these empirical t-matrices off-shell 
without using a realistic potential model. 

These LF on-shell t-matrices are actually energy de- 
pendent. In Fig. la, the t-matix involved is that for the 
initial dibaryon state of mass m* . We therefore use the 
empirical t- matrix at the N 2 energy of \fs = m*. 

Now the d* — > NN decay in free space gives two nucle- 
ons each of momentum p max - The resulting decay width 
is represented by a solid circle in Fig. 3. In Fig. la, 
on the other hand, we need decay widths at the same 
\fs = m* but with off-shell nucleon momentum p < p max . 
They can also be calculated for the same t-matrix input 
from the decay width formula of (!]] because this formula 
works for any nucleon momentum. The result is shown 
as a dashed curve in Fig. 3, one for each y/s whose value 
(in MeV) is also given near the curve. These t-matrices 
are said to be calculated at the (initial-state) dibaryon 
energy (or DBE). 

In contrast, the recoil interaction in Figs, lb or lc 
takes place after pion emission, i.e. at lower energies. 
It involves an off-shell t-matrix defined at the energy of 
the final N 2 state, as in the analgous processes in pp 
bremstrahlung ]T5|-|l7|] . In the present rough estimate, it 
is convenient to interpret it as an off-shell BB — > NN 
decay process at an intermediate i? 2 -state energy y/s < 
m* . This means that one should use points from the 
dashed curves in Fig. 3 at smaller cm. energies. The 
figure shows that such decay widths are smaller because 
the empirical i-matrix is found to decrease rapidly as y/s 
decreases. 

Such a strong energy dependence originates from the 
strong energy dependence of certain empirical NN phase 
parameters such as the mixing parameter ei in the 
N 2 ( 3 Si — 3 Di) states. The LF t-matrices we use were 
constructed in [pi from the energy-dependent phase so- 
lution SP84 of Arndt and collaborators |Q. It is well- 
known that the energy-dependent mixing parameter ei 
from SP84 already varies more smoothly and more slowly 
with energy than the so-called single-energy solutions. 
Nevertheless, it still has an energy dependence much 
stronger than that of any common NN potential such as 
the FB potential. This situation can be seen in Fig. 15-7 
of ||. In other words, the empirical t-matrices and the 
FB potential do not represent exactly the same dynam- 
ical input when it comes to this tricky isovector tensor- 
force part of the NN interaction. 

Given this complication, and the other problem of 
missing off-shell dynamics, we decide to use a much 



cruder approach for these empirical t-matrices. Instead 
of using them at different effective B 2 energies dictated 
by each diagram, we shall use two extreme prescriptions 
meant to provide rough upper and lower bounds. 

The first prescription is to use t-matrices of the same 
NN energy yfs = m* in all decay diagrams, i.e. using 
only points lying on one of the dashed curves in Fig. 3. 
These t-matrices are thus calculated at the initial-state 
dibaryon energy, a prescription referred to below as the 
LF-DBE model. One can see from Fig. 3 that the re- 
sulting decay widths are larger than those calculated at 
variable energies. The DBE results for T^nn can there- 
fore be expected to be rough upper bounds. Fig. 3 shows 
that these DBE results are likely to be close to those cal- 
culated with the FB dynamics. 

The second prescription uses t-matrices calculated at 
the lowest possible energies in order to generate a rough 
lower bound. These energies will be taken to be NN 
energies (NNE) after the recoil, giving a prescription to 
be called the LF-NNE model. The resulting r 2a decay 
widths arc the free-space widths reported in |Q and rep- 
resented by solid circles in Fig. 3. In other words, the 
solid circles are now treated as a function of the nucleon 
momentum p2 and used in all diagrams. We can see from 
Fig. 3 that the T^nn that results will be much smaller 
than those for the LF-DBE or FB models. Since the 
LF-DBE model does not contain the effects of decreased 
effective B 2 energies for the recoil after pion emission 
while the LF-NNE model uses NN energies that are far 
too low, the actual LF results can be expected to be 
somewhere between these two LF models. 

Although we have not calculated the full off-shell f- 
matrix from the FB potential, it is nevertheless useful to 
point out some expected re-scattering effects in it. Wc 
shall concentrate on r 2a , where the A 2 ( 7 S 3 ) — » N 2 ( 3 D 3 ) 
matrix element of the isovector tensor force has the same 
sign as the N 2 ( 3 Si) — > N 2 ( 3 Di) matrix element, accord- 



ing to Eqs. (A3) and (A4). The long-range part of this 
isovector tensor force is from pion exchange and is attrac- 
tive in the N 2 ( 3 Si— 3 Di) states. Hence re-scatterings can 
be expected to enhance its contribution, i.e. to increase 
the width in Fig. 3 for small momenta. On the other 
hand, the short-range part of the tensor force is domi- 
nated by rho exchange, and is repulsive. Re-scatterings 
will then reduce its value, leading to a decreased width in 
Fig. 3 at large momenta. These considerations suggest 
that with re-scatterings, the width from the FB model 
is likely to move towards that for the LF-DBE model. 
Even though we do not know the quantitative extent of 
the change, we do expect the actual results to be much 
closer to the LF-DBE or FB values than the LF-NNE 
values. This is especially true if the mixing parameter 
ei given by the FB potential is more realistic than the 
empirical values obtained from phase-shift anaylsis. 

To summarize, the three models shown in Fig. 3 seem 
to present an interesting range of dynamics for our rough 
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estimate of the nNN decay width. 

Let us now return to our rough estimate of T a using 
the FB model for T 2a . We average the dimensionless "re- 
duced" width T 2a /p 2 by itself over a rectangular distribu- 
tion of p 2 (for approximately nonrelativistic kinematics) 
to get the rough estimate 



T2a 
P2 



1 



I'm 



^2ma 
0.010 (0.022) 



d(j? 



(39) 



for m* = 2100 (2200) MeV. 

The energy denominator D a is next factored out of the 
integral by using the median energy (E±) of E\ given in 
Eq. ((2^). After replacing the factor E\ in the integrand 
itself by its median value, we are finally left with an in- 
tegration over the E%E 3 (or Nn) Dalitz area ju| 



A= I dE 2 dE 3 . 

'Dalitz 



(40) 



This area turns out to be 1390 (6800) MeV 2 . Hence 



25_ A /T 2a \ {Eu 



- 36 



12tt \p 2 / (D a )i \ p 3 
0.04 (0.7) MeV, 



(41) 



for the FB potential. As expected, the three-body de- 
cay width increases sharply with increasing to*. It is 
also rather small, even before the destructive interference 
with the other amplitudes is included. 

It is now a simple matter to return to Eq. (Pq) and 
perform an honest integration over the Dalitz area with- 
out making the separate averages decribed previously in 
this section. In the energy denominators, we use the ac- 
tual external energies E(pi) and E(p 2 ) of each point in 
phase space, the approximate internal energies for Ei(p) 
shown in Eqs. (^0|) and (|3l|), and the approximation 



Ei(p-p 3 ) 



<p z > 



Pi 



(42) 



that is the leading term of a Legendre expansion in p.p3. 
The calculated results from Fig. la alone are shown in 
Fig. 4 as a light solid curve for the FB model, as open 
diamonds for the LF-DBE model, and as a light dashed 
curve the LF-NNE model. The results for Figs, la-c are 
shown as a heavy dotted curve for the FB model, solid 
squares for the LF-DBE model, and a dot-dashed curve 
for the LF-NNE model. 

We see that the results from all three models are in 
better agreement the higher the d* mass. The FB and 
LF-DBE results are quite close to each other over the 
entire mass range, with the LF-DBE widths a little larger 
close to the threshold, and a little smaller as to* increases 
beyond 2100 MeV. The LF-NNE widths are about an 
order of magnitude smaller than those from the other two 



models below 2100 MeV. Hence the FB widths seem to 
be quite reasonable, being comfortably within our rough 
upper and lower bounds at low d* masses. 

Figure 4 also shows clearly how the destructive interef- 
erence among Figs, la-c reduces their resultant decay 
width significantly, as we have already noted previously. 



V. RETARDATION AND PAIR 
CONTRIBUTIONS 

The remaining six diagrams shown in Fig. 1 come in 
two groups: Figs, ld-f describe the retardation contribu- 
tions due to pion emissions from a baryon when a virtual 
meson is in the air, while Fig. lg-i describe pion emission 
via the creation of a BB pair. 

The propagator for the retardation diagram a con- 
tains the energy-denominator factors Dq 1 D~ 1 , where the 
initial-state energy denominator 



D Q =m* -Lo-E(p 2 )-E A (p) 



(43) 



is the same in all three diagrams. Note that the sum of 
the two intermediate-state baryon energies that appears 
is only about 2210 (2240) MeV for to* = 2100 (2200) 
MeV. That is, this energy sum increases very slowly be- 
yond the sum of their rest energies (at 2170 MeV) as to* 
increases. Hence for to* as 2200 MeV, Do is close to — u>, 
where uj — y/ to 2 ^ + q 2 is the meson energy. 

The final-state energy denominator D a depends some- 
what on the diagram a with: 



D d = E(p 2 )-Lu- E A {p), 

D e = E{ Pl )-u-E{-p-p 3 ), 



D 



I 



E(pi 



E A (-p-p 3 ) 



(44) 



All angle-dependent energies will eventually be approxi- 
mated by angle-averaged energies. 

We turn next to the integration over one of the inter- 
nal momenta, say the meson momentum q. (The other 
internal momentum p is then fixed by momentum con- 
servation.) This integration works out the same way it 
does for the two-body transition amplitude T of Fig. la- 
c fO]. The major difference is that the meson propagator 
u> for Figs, la-c is now replaced by the propagator 
l/ujD()D a for a retardation diagram a. The appearance 
of an additional w-dependent energy denominator has the 
consequence that the longer-range 7r-exchange contribu- 
tion becomes relatively more important than in Figs, la- 
c. 

This change of propagator causes some complication in 
the calculation. Since the additional energy denominator 
depends on the meson mass, it can be included only if 
the 7r- and p-exchange contributions can be separated. 
The transition amplitude T from Figs, la-c is of course 
separable into the form: 
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(45) 



when the virtual-meson exchange is treated literally in 
the Born approximation implied by the figures. By the 
same argument, the retardation diagrams ld-f taken lit- 
erally as Born diagrams give rise to transition amplitudes 
1Z of the same separable form 



K = 1Z 7T +1Z P . 



(46) 



The situation is more complicated when re-scattering 
effects are also included, for example, by using the LF 
t-matrices. While the amplitude T can be calculated 
without the separation of terms shown in Eq. ((45[), the 
calculation of 1Z is not possible without it. Fortunately, 
each LF i-matrix is made up of a sum of three or four 
terms with different effective meson masses, although the 
function used differs somewhat from that appearing in 
one-boson exchange potentials. Hence the calculation of 
1Z is still possible if one interprets these terms as the 
separate contributions from different mesons or groups 
of mesons. 

This prescription is admittedly far from ideal, but even 
a theoretically calculated i-matrix may not be separable 
into terms describing single-meson exchanges. The rea- 
son is simply that re-scattering means the inclusion of 
terms with two or more exchanged mesons. Their cor- 
rect treatment can only be based on Feynman diagrams. 
In spite of this limitation, a calculation using the empir- 
ical i-matrix should still be very informative. The actual 
formula used to calculate 1Z for LF i-matrices will be 
given in the Appendix. 

One final problem in using the LF t-matrices has to 
be mentioned. The longest-range part of their tensor t- 
matrices is not one with the pion mass, but one with twice 
the mass. The authors found the constraint to match the 
one-pion range "too restrictive" though they were 
successful in constraining the longest-range term of the 
real part of the central t-matrices to the one-pion value. 
This range problem is probably not serious in generating 
the amplitude but possibly more serious in the ampli- 
tudes 1Z for retardation diagrams. This is because these 
diagrams contain one more energy denominator depen- 
dent on the meson energy u. If the meson mass is too 
large in these energy denominators, the resulting decay 
amplitudes might be too small. It is hard to tell if this is 
indeed the case for a numerically fitted function as there 
might be compensating terms in the fit, but we shall 
watch out for this possibility in the calculated results. 

We turn next to the pair diagrams of Figs. lg- 
i. Since we need a rank-3 spin operator for the decay 
d*[A 2 ( 7 S , 3 )] — > irNN, the pion-emission vertex has to be 
spin-dependent, while the virtual-meson exchange from 
or to the pair must involve a rank-2 spin operator. The 
spin-dependent part of the pion emission vertex with the 
creation or annihilation of a BB pair turns out to be of 
order p/m (ratio of baryon momentum to baryon mass) 



relative to the spin-dependent pion emission vertex from 
a baryon. The rank-2 spin part of the B 2 B — > B inter- 
action turns out to be also of order p/m relative to the 
rank-2 spin operator in the nonrelativistic BB — > BB in- 
teraction for both 7T and p exchanges. Furthermore, the 
additional energy denominator for these pair diagrams 
is larger in absolute value than those in Figs, la-f by 
roughly twice the baryon mass. The consequence is that 
these pair contributions to the decay amplitude can be 
expected to be roughly two orders of magnitude smaller 
than those for the other diagrams of Fig. 1 in the mass 
region of interest in this paper. They can therefore be 
neglected. 

Returning now to the retardation amplitudes 7Z, we 
need to add them to Eq. (|2^) and perform spin averages 
in the same way as for Figs. la-c. Again, it will be 
convenient to express all operator matrix elements for 
the retardation diagrams in terms of that appearing in 
Fig. Id for which the transition amplitude is 7Z 2d , where 
the subscript 2 refers to the "spectator" nucleon. The 
final expression for the irNN width for Figs, la-f turns 
out to have the same structure as Eq. (|25|): 



af 



where 



25 

127T 



\P2P3 
2 



B(a,d) + S p [ -B(b,e) + -B(c,f) 



(47) 



B{a,d)^1Z 2d {d;p 2 ) + 
B(b,e)^lZ 2d (e;p* 2 ) + 
B(c,f)tt1l 2d (f; P * 2 ) + 



F 2a (p 2 ) 

D a ' 

D b : 

^2a(p* 2 ) 



(48) 



are amplitudes for pion emission from a baryon. The di- 
agram labels a — f are kept in the expressions as they 
define the energy denominators to be used in the ampli- 
tudes. 

To simplify the calculation of the retardation ampli- 
tudes, we use in the energy denominators D e and D/ the 
approximation 



E t (-p - p 3 ) - E{-p 2 - p 3 ) » 
Ei(-p) - E(-p 2 ) 



(49) 



where i = A or N (i.e., blank). The mean external nu- 
cleon energy defined in Eq. (E8f) is used for E(±p 2 ), and 
the r.m.s. momentum (p 2 ) 1 ^of the initial state is used 
in the internal energies Ei(p). 

The resulting irNN decay width T a f is shown as a solid 
curve for the FB model, and as solid circles for the LF- 
DBE model. For comparison, the partial width T a for 
Fig. la alone is also shown, as a light solid curve for the 
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FB model, and open diamonds for the LF-DBE model. 
All these widths have been calculated for a d* radius of 
0.7 fm. The width T af for a d* radius of r* = 0.9 (0.5) 
fm is also shown in the figure as a long-dashed (dashed) 
curve. 

The main feature seen in Fig. 5 is that the retardation 
diagrams have overcome the reduction caused by the de- 
structive interference among Figs, la-c to give a more 
normal result quite close to r a for Fig. la alone. It is 
interesting that T a / is greater than r a for the FB model, 
but smaller for the LF-DBE model. The difference has a 
simple explanation. 

First of all, we should recall that the decay at the 
BB — > NN step is controlled by a competition between 
the long-range attractive tensor force from 7r-exchange 
and the shorter-range repulsive tensor force from p- 
exchange |jj . This cancellation makes the width increase 
less rapidly with increasing dibaryon mass. It also ex- 
plains why the width is smaller for a dibaryon of smaller 
size, as seen in Fig. 5. 

Now the retardation diagrams contain an extra energy 
denominator that depends on the meson mass, a feature 
that favors the contribution from the virtual meson of 
smaller mass. Consequently, the attractive 7r-exchange 
contribution is re-inforced in the FB model. In the 
LF models, the one-pion term is absent in the isovec- 
tor tensor force. The terms that remain correspond to 
larger meson masses giving larger energy denominator, 
and hence smaller retardation contributions. 

If this picture is correct, the FB result would appear 
to be more sensible. However, a careful study of both 
theoretical i-matrices calculated from potential models 
and empirical i-matrices constructed from phase shifts 
will be needed to confirm this interpretation. We have 
not calculated the results for the TF-NNE model, but 
they are likely to be much smaller. 

Our general conclusion is that while the energy depen- 
dence of the i-matrix is not insignificant, the best esti- 
mate we have now is based on the FB model. According 
to Fig. 5, it gives a value of r a / for Figs, la-f that is 
about 0.006 (0.07, 0.5) MeV at to* = 2065 (2100, 2150) 
MeV. 

VI. I FORBIDDENNESS AND THRESHOLD 
BEHAVIORS 

It is now useful to discuss the threshold behavior of 
the decay width. We shall begin with the Dalitz area 
A(m*) defined in Eq. (|o|). In the NR limit, it varies 
with momenta roughly as (p2maxP3max) 2 - In reality, the 
dependence is only a little weaker than this, as the dotted 
curve for the ratio 

CA(m*) = A(TO*)/(p 2 maxP3ma X ) 2 (50) 

of A to its threshold behavior shows in Fig. 6. 



Using the NR dependence for the Dalitz area, we find 
a threshold behavior for the 3-body decay of 

T(m*) w const pf^™- ( 51 ) 

This dependence could be approximated by just Pamaxi 
where £ = £^+£n is the degree of forbiddenness of the de- 
cay, but we shall not make this additional approximation 
here. 

For d* — > itNN with P-wave pion emission, we have 
£ n = 1, In = 2, and a decay that is third forbidden. To 
see how dominant this threshold power law is, we show 
in Fig. 6 the ratios 

C a (m*) = r Q (TO*)/(^ max p^ max ), (52) 

where a = a and af, as functions of the d* mass for both 
the partial width r a from Fig. la only (dashed curve) 
and the width r a / from Figs, la-f (solid curve). Both 
widths are calculated with the FB model of interactions. 

These curves show that the threshold behavior is fur- 
ther modified by momentum dependences due to dynam- 
ics and form factors. In fact, both decay widths behave 
roughly as P3 max P2max over the first half of the mass 
range shown in the figure, a small reduction in the mo- 
mentum power being used to simulate the significant mo- 
mentum dependences shown in the figure. 

These results are relatively simple because the initial 
state is a relative S-state. The situation is more compli- 
cated for a relative D-state. The solid spherical harmonic 
3^2m(p) that appears in the initial-state wave function 
would have given rise to an additional factor of 

P 2 = (p 2 + q) 2 , (53) 

where q is the momentum of the virtual meson involved 
in the BB recoil. Since this internal momentum is even- 
tually integrated over, what survives of the q 2 term has 
the same threshold behavior as those from the S-state 
component of d* . The remaining terms in the equation 
have additional powers of the external momentum p2 left 
over, leading to amplitudes with higher powers of p2- 

We have shown in Sec. II that the d* decay by S-wave 
pion emission is also third forbidden. This S-wave emis- 
sion comes from the second, or S-wave, term in the basic 
irqq vertex calculated from the usual irqq Lagrangian with 
either pseudovector or pseudoscalar coupling: 

V wqq = const[-<r.(p - p') + (w/4m)cr.(p' + p)] (54) 

in the NR and w << to limit Here p (p') is the quark 
momentum before (after) pion emission, the momentum 
of the emitted pion being P3 = p p', and w = E — E' 
is the baryon energy transfer. One can show that for 
d* decay, the use of the second term, the S-wave irqq 
vertex, instead of the first, or P-wave, vertex results in 
the substitution of one power of the pion momentum p% 
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by one power of the NN relative momentum |pi — P2 1 
in the final state. This substitution leaves the threshold 
power law unchanged. 

However, it is not true that the threshold power law 
is always determined completely by the degree of forbid- 
denness of the decay. Occasionally the expected thresh- 
old term is absent. Then the surviving leading terms 
are higher-order Taylor terms with a momentum power 
greater than expected. This situation will be described 
as being "unfavored" because an increase in the power 
implies a reduction in the decay width near threshold. 
So the general rule for the threshold power law is that 
the power is not smaller than that determined by the l- 
forbiddenness of the decay, but it could be bigger if the 
decay is unfavored. 

To support the assertion that some decay amplitudes 
are unfavored, I now describe an example which is in- 
teresting in its own right. It is the irNN decay of 
that most promising of low-mass dibaryons, namely d! 
(J n = 0~,T = or 2), which was first proposed for a 
dibaryon interpretation of the structures seen in the pi- 
onic double charge exchange reaction nn{n + 1 'K~)pp on 
nuclear neutrons 19-2j|], 

The decay equations for S- and P-wave pion emissions 
are, respectively 

d'(0-,T = 0,2)^iV 2 ( 1 S*o)+7r(£ 7r =0). (55) 

d'(0-,T = 0,2)^iV 2 ( 3 P 1 )+ 7 r(^ = 1); (56) 

This shows that S-wave pion emission is I = allowed, 
while P-wave pion emission is £ = 2 forbidden. This 
classification suggests that S-wave pion emission would 
normally dominate the decay near threshold. 

A more detailed analysis shows that this is not always 
true. To see this, let us start with the simple models of 
d! used in p2fl : a AN bound state in a relative P-state 
for the T — 2 dibaryon, and an N*N relative-S bound 
state involving a P-state N* for the T — dibaryon. 

In the AN model of d'(T = 2), the off-shell A can de- 
cay into a nucleon while the spectator nucleon can come 
out without further interaction. The pion emission ver- 
tex can be written in one of the following two forms when 
specialized to the three-body final state 



NA 



const[er.p 3 + (iy/2m)<T.pi 2 ] , 



VnNA — COnSt 

where 

P3 = Ptt 



/ w \ w 

V Am) 2m ^ 2 



Pl2 = ^(Pl - P2) = -P2 - ^P3 



(57) 



(58) 



(59) 



in the dibaryon rest frame, and where <x now stands for 
a transition spin operator M. We shall call these ex- 
pressions the Jacobi form and the single-particle form, 
respectively. 



The spectator nucleon (of momentum p 2 ) is a real 
spectator when final-state interactions are neglected. Its 
P-state wave function is 



-01m(P2) = /l(P2)3 ; im(P2), 



(60) 



where the radial function f\ is finite at p 2 = 0. 

Using the single-particle form for the pion emission 
vertex, we find that the first, or P-wave, vertex gives a 
decay amplitude proportional to 3^im(P2)o'-P3; while the 
S-wave vertex gives one proportional to 3 /, im(P2)f-P2- 
One can recognize the presence of a P-wave pion in the 
first decay amplitude, and an S-wave pion in the second 
amplitude. However, both amplitudes depend on the ex- 
ternal hadron momenta to the same second power. 

It is easy to see why the naively-expected threshold 
term for S-wave pion emission is missing. The S-wave 
vertex is actually a three-body interaction that is P-wave 
in the spectator momentum p 2 . Multiplication into the 
initial P-state wave function gives an extra factor of p\. 
In other words, the normal S-wave result at threshold, as- 
sumed to be momentum- independent, is simply absent. 
What survives is a higher-order term. We shall call at- 
tention to this abnormal situation by calling the decay 
unfavored. 

An amusing complication in this decay of the AN 
model d'(T = 2) is worth pointing out. In terms of the 
Jacobi momentum pi 2 of the final A*" 2 state, the initial 
P-state becomes a mixture of P- and S-states: 



3>lm(P2) = 3W-P12 - ^Pa) 

= -3WPl2) - ^lmfes) 



(61) 



This is caused by the recoil of the NN system on pion 
emission. The second term on the right-hand side de- 
scribes the process where the pion carries away the or- 
bital angular momentum leaving the ./V 2 pair in a relative 
S-state. 

The Jacobi description gives four amplitudes: Two are 
similar to those already described: a term proportional 
to 3^im(Pi2) (T -P3 for P-wave pion emission, and a term 
proportional to 3^i m (Pi2) cr -Pi2 f° r S-wave pion emission. 
However, two additional terms appear: 3^im(P3)f-Pi2 
describing P-wave emission from the S-wave vertex, and 
finally 3^i m (P3)"'-P3 describing S- and D-wave pion emis- 
sions from the P-wave vertex. These unexpected terms 
have their orbital angular momentum changed by baryon 
recoil on pion emission. The D-wave term comes with an 
AA( 1 D 2 ) pair in the final state. It is £ = 4 forbidden, 
and is therefore quite unimportant. These recoil-induced 
emissions are present in the single-particle form as well, 
but one has to know where to find them. 

One can see from Eq. (|57j) that the recoil-induced S- 
wave pion emission amplitude is not less important than 
the direct S-wave pion emission amplitude from the S- 
wave vertex, because unlike the latter it is not reduced 
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by the additional factor w/2m. However, it generates the 
same unfavored S-wave threshold behavior of a quadratic 
dependence on the external momenta and not the nor- 
mal behavior independent of external momenta. So none 
of the S-wave decay amplitudes in our simple model of 
d'(T = 2) is normal, and all are unfavored by two powers 
of the external momenta. 

The above analysis is helpful because it tells us how 
to generate S-wave emission amplitudes that are normal 
and therefore dominant near threshold: The unwanted 
momentum can be prevented from appearing as an ex- 
ternal momentum in an S-wave emission amplitude if it 
can be changed into an internal momentum that is even- 
tually integrated over. This means that the P-wave ex- 
citation of the initial state should not be in the relative 
BB coordinate, but in a quark coordinate internal to a 
baryon. In other words, amplitudes with normal S-wave 
threshold behavior can only come from those components 
of d'{T — 2) containing a P-wave excited A* or N* . 

The contributing term can readily be isolated by us- 
ing a shell-model configuration in which a P-state quark 
in B* becomes an S-state quark in a nucleon after pion 
emission. The angular integrations over the initial-state 
quark momentum p has the form 



rf 2 p3 7 im(p)K gg exp[-&(p-p 3 ) 2 ] , 



(62) 



where the solid spherical harmonic comes from the ini- 
tial P-state and the Gaussian comes from the final S-state 
wave function. Using the single-particle form Eq. (|5|) of 
the ttNN vertex with p2 = — p, we find that the P-wave 
vertex term contributes an S-wave pion emission ampli- 
tude that is unfavored by a momentum factor p|. It is the 
S-wave vertex term that gives rise to the sought-for am- 
plitude that is independent of external momenta. Such 
an amplitude is more favorable than P-wave emissions 
near threshold by two degrees of forbiddenness. It can 
be expected to dominate the decay width near threshold. 

It is true, however, that the S-wave vertex is weaker 
than the P-wave vertex by a factor w/2m, and, depend- 
ing on the dynamical model for d! , the B*B components 
of d'(T = 2) might themselves be weak. So the quanti- 
tative impact of this normal S-wave decay amplitude at 
the dibaryon mass of 2065 MeV is model-dependent. 

Exactly the same situation holds for d'(T — 0) when 
treated as an S-wave N*N bound state, where N* is a 
P-state baryon of isospin 1/2. To the extent that this 
component might actually be a significant if not major 
part of d'(T = 0), the decay width of this dibaryon can 
be expected to be actually larger than that for d'(T = 2) 
near threshold. 

The qualitative view on the decay of d! implied by 
our forbiddenness classification and threshold behaviors 
both agrees with and differs from those obtained in recent 
studies of the d! decay width in interesting ways. For the 
same models of d', we agree with Garcilazo p2| that the 



T = decay width should be larger than the T = 2 
width, but both our widths are likely to be much smaller 
than Garcilazo's results. This is because for d'(T = 0), 
the normal allowed S-wave decay comes only from the 
weak S-wave nqq vertex, while for d'(T — 2) without any 
B* component, the I = 2 forbidden decay width should 
be relatively small near threshold. 



We agree with Obukhovsky et al. 23 in finding that 
the P-wave ttNN vertex causes a recoil-induced decay 
into S-wave pion and NN^So) final state, and that the 
resulting decay amplitude has an abnormal quadratic de- 
pendence on the pion momentum. However, another de- 
cay amplitude of comparable importance, namely P-wave 
pion emission caused by the P-wave nNN vertex and 
leading to the NN( 3 Pi) final state has not been included 
in their calculation. 

We agree with Samsonov and Schepkin p4j that an- 
other problem with the calculation of [^3| is that the 
S-wave irNN vertex has been left out. Its inclusion is 
important as it gives rise to the dominant allowed S-wave 
emission near threshold. Much is unclear, however. The 
energy transfer w appearing in it is sometimes positive 
and sometime negative. It is not clear what it averages 
out to be. It is also not clear how much of the S-wave 
dominance remains by the time we get to the dibaryon 
mass of 2065 MeV. 



VII. DISCUSSION AND CONCLUSIONS 

Just as in nuclear (3 decay, a general three-body decay 
near threshold depends sensitively on the orbital angular 
momenta of the decay particles in the final state. A clas- 
sification of the decay widths based on different degrees 
of ^-forbiddenness then becomes useful. Much diversity 
exists within this broad classification, primarily because 
some decays are weaker than normal (or are unfavored) 
due to the presence of additional powers of external mo- 
menta in their decay amplitudes. 

According to this classification, the decays 

d*(J* = 3+,T = 0) -> tt(4 = 0)NN( 3 F 3 ) 

and7r(4 = l)NN{ 1 D 2 ) (63) 

of the dibaryon d* are both normal "third- forbidden" de- 
cays. Even among such normal decays, considerable di- 
versity remains in the decay widths due to differences in 
wave functions and in dynamics. 

In this paper, we have used an angle-average approxi- 
mation to make quick estimates of the irNN decay width 
of a d* described as a simple A 2 bound state. We find a 
leading-order decay width of only 70 keV at the d* mass 
m* of 2100 MeV. The calculated width decreases rapidly 
as m* decreases towards the threshold at 2020 MeV, be- 
ing only about 6 keV at m* = 2065 MeV. 
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The accuracy of these angle-averaged results must be 
confirmed in the future by more detailed calculation. 
Such a calculation is conceptually simple, but a little 
tedious in execution as it involves additional angle inte- 
grations. The uncertainties caused by angle averaging are 
likely to be much smaller than the uncertainties caused 
by uncertainties in the wave functions and in the strong- 
interaction dynamics. 

Although the d* model used here is crude, the calcu- 
lated results can be used to estimate the decay widths 
for other models of the wave function. For example, the 
d* wavefunction is much more complicated in the quark- 
delocalization and color-screening model of It has 

been estimated |fflj that quark derealization would de- 
crease the present result by a factor of 0.4, while only 1/5 
of the six-quark wavefunction is in the A 2 configuration. 
The remainining 4/5 of the d* state is in hidden-color 
configurations which will contribute much less, perhaps 
only 1/5 or less of the contribution of A 2 . 

The d* — + NN width used here might also have been 
overestimated by the meson-exchange model used in [0. 
This would be the case if the interior of the d* dibaryon 
is a perturbative vacuum where the exchanged mesons 
do not exist as such and where the dynamics is much 
weaker. Assuming that the hidden-color contribution is 
cancelled by the reduction caused by the perturbative 
vacuum, we are left with a total reduction, for the quark- 
delocalization and color-screening model of d* , of about 
an order of magnitude from the results reported here for 
our simple model of d* as A 2 . 

There are other wave- function uncertainties such as the 
presence of the A 2 ( 7 Da) component, which could be sig- 
nificant or negligible depending on the dynamics assumed 
to operate at short distances. 

We have relied on the familiar meson-exchange model 
of nuclear forces to generate the dynamics needed in the 
calculation. Unfortunately, the calculated decay width 
involves the short-range part of the isovector tensor force 
where there is much cancellation between the tt- and p- 
exchange contributions. Our knowledge of this particu- 
lar combination of isovector tensor forces, and of short- 
distance dynamics in general, seems to be relatively poor. 

In spite of these uncertainties, our results seem to sug- 
gest that the nNN decay width of d* might be too small 
at low d* masses to make it easy to use the few decay 
pions for particle identification. 



APPENDIX A: SPIN AVERAGES AND 
TWO-BODY DECAY AMPLITUDES 

The three-body decay d* — > ttNN involves two groups 
of two-body reduced matrix elements (RME) - one for 
7r emission from a baryon and one for the two-baryon 
transition BB — > NN, where a baryon could be A or 
N. When brought down to the quark level, the RME's 
needed for tt emission from different baryons are related 
by the nonrelativistic (NR) quark model as follows: 

(AUffmlltf) = jp/2, 

(A||triTi||A) = y, (N\\<r 1 T 1 \\N) = ^. (Al) 

where quark 1 is in the baryon. 

The operator responsible for the BB — > NN transi- 
tions needed in Figs, la-c is 



BB 



-(<Tl X CT 4 ) (2) (Ti.T4) 



(A2) 



where quarks 1 and 4 are in two different baryons. The 
following RME's are needed in our calculations: 



(A 2 ,ST = 30||T B b||A 2 ,10) 
(JVA, ST = 21||T BB ||A 2 ,01) 
(A 2 ,ST = 21||T BB ||A 2 ,01) 



16 7 
TV 2' 
40 



27 



V5, 



40 
27 



y/2 



(A3) 



For comparison, we also give the RME between N 2 
states: 

50 

(A 2 ,ST=10||T BB ||A 2 ,10) = -\/5. (A4) 

The spin averages of the product of operator matrix el- 
ements appearing in the contributions to the decay width 
shown in Eq. (|l^) are relatively easy to calculate directly 
for the d* S-state. For Fig. la alone, we have 



2S + 1 ' ' 

M,M T fi,m,n 

25, 



(i\(cri x cr 4 )( 2 )(i-i.T 4 )cr m T„|/)| 2 



= ^(A 2 , ST = 30||T BB || A 2 , 10)(A||<titi|| A) , (A5) 
where 
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\i) = \d*;S = 3M,T= 00), 
|/) = |(7rAW)00,LM T ) • 



(A6) 



The quark label has been left out of the quark operators 
VmTn responsible for the pion emission. We have also 
chosen to express this and other spin averages in terms of 
the RME's for A 2 — > NN and A — > irN for convenience 
in mutual comparisons. The spin averages for the other 
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diagrams and for interference terms can be calculated in 
the same way. 

The pion momentum is the same for all Feynman di- 
agrams, but the dependence on the nucleon momenta is 
more complicated. We shall assume that the BB — > NN 
interaction should be evaluated in the cm. frame of the 
baryons in Figs, lb, c, e and f, and not that of d* . This 
means, for example, that the nucleon momentum appear- 
ing in the corresponding decay amplitudes for Figs, lb 
and c, denoted T 2a in the text, is different from that for 
Fig. la, as explained more fully in Sect. II. This refine- 
ment is included in the final expression p7|), which also 
gives the results of all spin averages. 

Two two-body decay widths (or amplitudes) appear in 
this equation. The one for the decay A — > irN is obtained 
from the free-space decay width Ta based on a P-wave 
irqq vertex for pion emission from quarks. This means 
that we have neglected a certain well-known S-wave ver- 
tex term which when expressed in baryon coordinates 
is proportional to the mean momentum of the emitting 
baryon. This S-wave emission vanishes in free space be- 
cause the appropriate frame to be used is the one where 
the mean momentum vanishes, i.e. the Brcit frame. It 
does not vanish in a many-baryon system where the mean 
baryon momentum is nonzero. We have argued in Sec. 
II, however, that for d* decay this S-wave vertex is likely 
to be much less important than the more familiar P-wave 
vertex. It may therefore be neglected in the qualitative 
estimates attempted in this paper. 

Two different dynamical inputs for the BB — > NN 
transition amplitudes T 2a are used in Eq. (ft7j). It de- 
pends on the nucleon momentum p\ in the B 2 rest frame 
of invariant mass m\ 2 : 




1 



7T/3* 



3/4 



where 



m 12 /4, K=p* 2 /P*, /?* 



(A7) 



(A8) 



r* is the d* radius, usually taken to be 0.7 fm, and I(k) 
is an integral jlj over the momentum q = [3*Q of the 
virtual meson responsible for the BB recoil: 



- K 2 /2 



-Q 2 /2 



[j 2 {inQ){nQf] 



t ivt (P*Q)QdQ- 



(A9) 



tivt(q) ~ v ivt (q) 



Vi(q), 



where 



Vi(q) = 4tt 



SjOj 1_ 

Am 2 q 2 + j 



A?-mr 2 

A? + q 2 



(A10) 



(All) 



The sum is taken over the two virtual mesons i = 7r, p, 
with signatures S„ = — 1, S p = 1 and strengths 

a„ = = 14.4, 

a P = ^( 1 + —) = 0.84(1 + 6.1) 2 , (A12) 
4tt V g v J 

respectively. The other parameters are those shown on 
p. 37 of @. 

For the Love-Franey i-matrix, we use Eq. (15c) of pa] 



tivt(q) = ^2u(q), 



where 



U(q) = 327T 



V?q*R 7 



[1 + {qRi 



,213 ■ 



(A13) 



(A14) 



This is a sum of 3-4 terms of different ranges Ri cor- 
responding to virtual mesons of different masses TOj = 
hc/Ri. The parameters are given in Table I of ||. 

In the approximation described in Sec. V, two of the 
retardation amplitudes 1Z involve the same function: 



(A15) 



Only two distinct functions are then needed. Both are 
generated from Eq. (A7), but by using a certain retar- 
dation analogs r ivt of t ivt . For the FB potential, they 
are 



nvt(a;q) = 



DoiD a i ' 



(A16) 



where the argument a = d (e) denotes Fig. Id (le), 
u>i = y/ m 2 + q 2 , and 



D 0l w m* - Ui - (E 2 ) - (Ea), 



D r 



(E 2 



{El 



(A17) 



The baryon energy Eb is Ea for a = d, and the nucleon 
energy E for a = e. The averages appearing here are 
those described in Sec. V. 



For the LF i-matrices, the Born term Vj, in Eq. ( |A1G| ) 
should be replaced by the f-matrix term tj. 



The function tivt (g) is the radial part of the isovector 
tensor force appearing in both meson-gg and meson- BB 
interactions. 

For the Full-Bonn (FB) potential Q in the Born ap- 
proximation, we use 



[1] C.W. Wong, Phys. Rev. C 57, 1962 (1998). 



14 



[2] T. Goldman, K. Maltman, G.T. Stephenson, Jr., K.E. 
Schmidt, and F. Wang, Phys. Rev. C 39, 1889 (1989); F. 
Wang, J.L. Ping, G.H. Wu, L.J. Teng, and T. Goldman, 
ibid 51, 3411 (1995). 

F. Wang, G.H. Wu, L.J. Teng, and T. Goldman, Phys. 
Rev. Lett. 69, 2901 (1992). 

T. Goldman, K. Maltman, G.T. Stephenson, Jr., J.L. 
Ping, and F. Wang, "Systematic Theoretical Search for 
Dibaryons in a Relativistic Model", Los Alamos Report 
LA-UR-95-2609. 

R. Machleidt, K. Holinde, and Ch. Elster, Phys. Rep., 
Phys. Lett. 149C, 1 (1987). 

M.A. Franey and W.G. Love, Phys. Rev. C 31, 488 
(1985). 

H. Sugarwara and F. von Hippel, Phys. Rev. 185, 2046 
(1969). 

T. Ericson and W. Weise, Pions and Nuclei (Clarendon 
Press, Oxford, 1988), pl20. 

T.-S. H. Lee and D.O. Riska, Phys. Rev. Lett. 70, 2237 
(1993). 

C. Hanhart, J. Haidenbauer, A. Reuber, C. Schutz, and 
J. Speth, Phys. Lett. B 358, 21 (1995). 
H.T. Coelho, T.K. Das, and M.R. Robilotta, Phys. Rev. 
C 28, 1812 (1983). 

Particle Data Group, Review of Particle Physics, Phys. 
Rev. D 54, 1 (1996). 

G. Kallen, Elementary Particle Physics (Addison- Wesley, 
Reading, MA, 1964), pp. 193-199. 

G. E. Brown and W. Weise, Phys. Rep., Phys. Lett. 22C, 
279 (1975). 

V.R. Brown, Phys. Rev. 177, 1498 (1968). 
V. Herrmann and K. Nakayama, Phys. Rev. C 46, 2199 
(1992); F. de Jong and K. Nakayama, Phys. Lett. B 385, 
33 (1996). 

J. A. Eden and M.F. Gari, Phys. Rev. C 53, 1102 (1996). 
R.A. Arndt, L.D. Roper, R.A. Bryan, R.B. Clark, B.J. 
VerWest, and P. Signell, Phys. Rev. D 28, 97 (1983); 
R.A. Arndt and L.D. Roper (unpublished); R.A Arndt, 
private communication via G. Love. 

R. Bilger et ai, Z. Phys. A343, 491 (1992); R. Bilger, 

H. A. Clement, and M.G. Schepkin, Phys. Rev. Lett. 71, 
42 (1993); 72, 2972 (1994). 
L.S. Vorobev et al, JETP Lett. 59, 77 (1994). 
W. Brodowski et ai, Z. Phys. A 355, 5 (1996). 

H. Garcilazo, J. Phys. G 23, 1101 (1997). 

I. T. Obukhovsky, K. Itonaga, Georg Wagner, A.J. Buch- 
mann, and Amand Faessler, Phys. Rev. C 56, 3295 
(1997). 



th/9712079 



nucl- 



A. Samsonov and M. Schepkin, Los Alamos Eprint 

M.A. Franey, Phys. Rev. C 24, 1073 



FIG. 3. The d* -> AW decay width parameter V 2a (p)/p 2 
as a function of the off-shell nucleon momentum p in the NN 
cm. frame for the Full Bonn (FB) and the Love-Franey (LF) 
models of dynamics. Each LF curve is labeled by the cm. 
energy of the NN system whose empirical free-space f-matrix 
is used in the calculation. The on-shell decay width at each 
energy is given by a solid circle in the figure. 



FIG. 4. The widths F a from Fig. la alone and F ac from 
Figs, la-c for the decay d* nNN for different dynamical 
inputs: the Full Bonn potential (FB), Love-Franey t- matrices 
calculated at the dibaryon energy (LF-DBE) and at the NN 
energy (LF-NNE). 



FIG. 5. Ratios of the Dalitz area A, the decay width r a 
from Fig. la, and the decay width r„/ from Figs, la-f to 
the expected threshold power law for the d* — > nNN decay 
as functions of the dibaryon mass m*. Each ratio is given in 
units of the value at m* = 2034 MeV. 



W.G. Love and 
(1981). 



Figure Captions 

FIG. 1. Some leading-order Feynman diagrams for the 
decay d* —* nNN with pion emission from a baryon. 



FIG. 2. Some leading-order Feynman diagrams for the 
decay d* — » nNN with pion emission from virtual mesons. 
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